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We obtain the following analytical formula which describes the dependence of the electric 
potential of a point-like charge on the distance away from it in the direction of an external 
Q_i' magnetic field B: 



r-| ' <&(z) = e/\z\ 1 - exp(-A/6m2 \z\) + exp(-A/(2/7r)e 3 _B + <om? e \z\) . 

The deviation from Coulomb's law becomes essential for B > 3i\B cr ja = 3itm 2 e /e 3 ~ 6 • 10 16 
G. In such superstrong fields, electrons are ultra-relativistic except those which occupy the 
CNj ; lowest Landau level (LLL) and which have the energy = ml +p 2 z . The energy spectrum on 
which LLL splits in the presence of the atomic nucleus is found analytically. For B > 3ixB cr / a 
it substantially differs from the one obtained without accounting for the modification of the 
atomic potential. 

1 Introduction 

^ ■ In the pioneering papers [Tj an exponential modification of Coulomb's law by superstrong 
magnetic fields B > m^/e 3 = 137B cr was discovered (B cr = m 2 e /e = 4.4 • 10 13 G = 4.4 • 10 9 T 
is the so-called critical or Schwinger magnetic field 0). It originates from the strong (linear) 
dependence of the photon polarization operator on the external B. The Coulomb potential 
gets modified at distances l/m e > r > 1/V e 3 B = an/e, where an is the Landau radius (the 
size of the ground state electron wave function in the direction transverse to the homogeneous 
magnetic field). In [1] the shape of the modified potential was determined numerically. 

As found long ago in papers [2] the photon polarization operator in a strong magnetic 
fields B > B cr is dominated by electron and positron states belonging to the lowest Landau 
level (LLL) and it factorizes into transverse and longitudinal (with respect to the direction 
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of B) parts. Its dependence on the longitudinal and time-like components of the momentum 
coincides with that in two-dimensional QED (with no B). In paper [3] a simple interpolating 
formula for the photon polarization operator in two-dimensional QED was found, which 
allows to find an analytical expression for the potential of a point-like electric charge in 
D = 2 QED. For light "electrons" propagating in the loop (m < g, where g is the electric 
charge in D = 2 QED), the electric potential gets screened at all distances larger than l/(2g). 

According to [3] the deviation of the interpolating formula from the exact expression for 
the polarization operator in D = 2 does not exceed 10%. In Section 2 we will show that for 
various values of the ratio g/m the accuracy of the analytical formula for the potential in 
D = 2 is always better than 4% for all distances z. 

For the realistic case D = 4 only the asymptotics of the electric potential at large (\z\ ^> 
l/m e ) and small (\z\ <C l/m e ) distances were found in [3]. In agreement with [1] the screening 
only takes place at short distances \z\ <C l/m e . In Section 3, using the interpolating formula 
for the polarization operator of [3], we will obtain an analytical expression for the dependence 
of the electric potential of a point-like charge on the distance to this charge in the longitudinal 
direction (parallel to the magnetic field B) and at zero transverse direction. This potential 
determines the energies of atomic electrons. The spectrum of electrons originating from LLL 
will be found in Section 4. 

In [1] the ground state energy of a hydrogen atom in a superstrong magnetic field was 
found in the shallow-well approximation. The detailed analytical study and comparison with 
numerical results for the atomic levels in an external magnetic field performed in paper [4] 
clearly demonstrates that the shallow-well approximation used originally for this problem 
in textbook [5] has a very poor accuracy (better to say no accuracy at all). An algebraic 
expression for the energies of the ground and excited states was obtained in |4] which repro- 
duces the results of numerical calculations with high accuracy. In Section 4 an analogous 
expression, but which now takes screening into account will be derived. It yields a value of 
the ground state energy in the limit of infinite B which is Eq = —1.7 keV, which strongly 
differs from the result of the shallow- well approximation [I]: Eq w = —4.0 keV. 

Excited states can be classified with respect to the change of sign of the z coordinate and 
divided into even and odd states. Odd states in the magnetic field B > ml/e 3 = B cr /a follow 
unperturbed Coulomb levels with very high accuracy: E odd = — me 4 /(2n 2 ), n — 1,2, ... and 
the screening does not alter this result. Concerning even states, at finite B they considerably 
deviate from Coulomb levels. These deviations decrease with increasing B like l/log(S), 
and, in the limit of infinite B, the even and odd states become degenerate j6]. The screening 
changes this result qualitatively and lifts this degeneracy: the energies of even states remain 
considerably higher than those of corresponding odd states, even in the limit B — > oo. 

The dependence of the electron mass on the external magnetic field is discussed in Ap- 
pendix. 

We summarize our results in the Conclusion. 

2 Coulomb potential modification in D = 2 QED 

The fermionic part of the D = 2 QED Lagrangian is L = ifj(id — gA — rnjip, where 70 = 
ci, 71 = 102 are Pauli matrices. Inserting the photon polarization operator II into the photon 
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propagator D leads to the following equations (see Fig. 1). The tree level Coulomb potential 
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where the expression for the photon polarization operator in D = 2 should be used. Instead 
of calculating the fermion loop we can take the expression for II obtained in the dimensional 
regularization method [7], substitute D = 2 in it and divide it by two, because in two 
dimensions the traces of 7-matrices are proportional to 2 instead of 4: 



u(k 2 ) = v 



\n(VT+t + Vi) - 1 



-Ag 2 P{t) 



(4) 



where t = —k 2 /4m 2 . 

This approach resembles the one used in the calculation of the Uehling-Serber corrections 
to the Coulomb potential in [8]. However we resum the whole geometric series instead of 
considering 1-loop correction. 



+ 
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Fig. 1. Modification of the Coulomb potential due to the dressing of the photon propagator. 



To obtain the electric potential &(z) of a charge at rest we should set k^ = (k ,k\\) = 
(0, fc||) (the notation k\\ is convenient in the D = 4 case) and, then, to make the Fourier 
transformation: 

e ik w z dkJ2TT 



$ z) 



Ang 



k 2 + Ag 2 P{k 2 /4m 2 ) 

—00 

Finally, the potential energy of the charges +g and —g is 

V{z) = -g<S>(z) . 



(5) 



(6) 



We did not succeed to analytically perform the integration in ([5]) with P(t) given by its 
exact expression (jlj); instead, an interpolating formula for P(t) was suggested in [3]. The 
asymptotics of P(t) are: 



such that the following interpolating formula 



P(t) 



2t 
3 + 2t 
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has a correct behavior at small and large t. Substituting (jHJ) in (JSJ and performing analyti- 
cally the integration yields [3] [§: 
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For heavy fermions (m 3> g) the potential is given by the tree level expression; radiative 
corrections are small: 



m ^> g 



In the case of light fermions (m <C 

m <C g 



p) we get a much more interesting behavior: 
-2tt(7 
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At short distances the potential has a Yukawa behavior, corresponding to the exchange 
of a massive vector particle. The fact that the photon remains massless for nonzero m 
follows from the behavior of the potential at large distances: a linear behavior corresponds 
to a massless exchange. However, the coupling constant differs from that at tree level by 
the small factor (3m 2 /2g 2 ): a screening of the tree level potential occurs at \z\ > l/(2g). 
Finally, for m = at all z the potential is given by the first line of ffTTj) : the photon gets a 
mass m 7 = 2g = 2^/a = e/y/n (where, and only in this formula, e is the coupling constant in 
the Heaviside units usually used in Quantum Field Theory) - a well-known result in D = 2 
QED with massless "electrons" (Schwinger model), where a nonzero mass for the gauge 
boson coexists with gauge invariance [9J[10J. 

We now study the accuracy of the analytical formula for the potential (Q compared with 
the numerical calculation of the integral (jSJ) in which the exact expression (Hj) for P(t) is 
used. In order to normalize the potential universally for the different values of the ratio 
g/m we choose an arbitrary constant up to which the potential is defined in such a way 
that &(z) = at z = 0. In Fig. 2 the exact (V(z)) and approximate (V(z)) potential 
energies are plotted together with their asymptotics for g = 0.5 and m — 0.1. The straight 
line which corresponds to the potential asymptotic at z — > oo crosses the vertical axis at 
V(0) w 7i-0(l - (3m/2#) 2 ) rs 1.42. 



3 We regularize the infrared divergence of the i integral making a subtraction at z = 0. 
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We see that the accuracy of the interpolating formula is very good. In Fig. 3 the relative 
difference (V — V") /V for the same values of g and m is shown. 
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Fig. 3. Relative difference of potential energies calculated with the exact and interpolating 
formulae for the polarization operator for g = 0.5, m = 0.1. 



The accuracy appears to be the worst at z ~ 1/m. We have checked that for various 
values of the ratio g/m, the relative accuracy of (Q is better than 4% at all z . 
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3 Screening of the Coulomb potential in D = 4 QED 
in a superstrong magnetic field 

In order to write an expression for the electric potential of a point-like charge we need a 
formula for the photon polarization operator in an external magnetic field B. This quantity 
was studied in many papers; a detailed review and references can be found in [IT]. In the 
calculation of the polarization operator, electron propagators in an external magnetic field 
are used. In the case of a strong field B ;§> B cr = m 2 /e the spacing of Landau levels is 
much larger than the electron mass. This is why the contribution of the electrons from LLL 
obeying the dispersion law = m 2 + p 2 dominates in the electron propagator. A linear 
dependence of the polarization operator on the external field B follows. Finally, the electric 
potential writes [HE]: 

*(*) = — 7 , (12) 

H + H + ^exp (-£) P (j| 

where k = (0, k x , k y , k z ), k\ = k 2 + k 2 , k z = k\\ and the magnetic field B is directed along 
the z axis. As stated in the introduction in (fT2l) a quasi two-dimensional formula for the 
photon polarization operator occurs. It was obtained in [2jj. Expression ( JT2l is valid for 
strong magnetic fields B 3> B cr = m 2 /e and for longitudinal momenta k 2 <C eB (the motion 
of the virtual electrons contributing to the photon polarization operator II occurs in one 
space and one time dimensions only for z larger than the Landau radius a# = 1/ V eB; for 
z < an the motion is four dimensional). 

Plugging in the Fourier transform of ( TT2|) the interpolating formula flH]) for P introduced 
in Section 2 for the electric potential of a point-like charge along the direction of magnetic 
field we get: 
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d(k\\\z\) , (13) 



where A is an ultraviolet cutoff in kj_ and as we will see later the values (kjjmax = A 2 < 2eB 
are essential. This condition allows us to replace the exponent in the denominator of the 
first line of f[T3"|) by one. This is why this exponent is absent in the second line of ffT3T) . 
Integrating by parts we get that, indeed, A 2 < k\ is important, and k\ eB for z 3> an. 
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So doing we obtain: 
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Performing the last integrals with the help of residues we finally obtain: 
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For magnetic fields B <C 37rmg/e 3 the potential is Coulomb up to small power suppressed 
terms: 
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in full agreement with the D = 2 case (see ([TO]) , where g 2 plays the role of e 3 B). In the 
opposite case of superstrong magnetic fields B 3> ?>nm 2 /e 3 we get: 
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which also closely resembles D = 2 case, ffTTj) . with the substitution e 3 B — > g 2 . 

The expression for the potential energy of two opposite charges separated by the distance 
z along the magnetic field (the transverse separation is zero) is: 



V(z) = -e*(z) = —. 



1-e" 



|»] + e - A /(2/ 7 r)e3S+6r, 



(18) 



In Fig. 4 the potential energy which follows from eq. f lTB]) for a magnetic field B = 10 17 
G is shown as well as its asymptotics at large distances: 
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Fig. 4. Modified Coulomb potential energy at B = 10 17 G (blue) and its long distance 
(green-pale) and short distance (red-dashed) asymptotics. 



In the following two figures (Fig. 5 and Fig. 6) the analytical formula (TIB]) for the 
potential energy is compared with the result of a numerical integration using the exact 
formula for P(k?J 'Am 2 e ) given by equation (jlj). The same value of the magnetic field B = 10 17 
G is used. We checked that for various values of B the agreement between analytical and 
numerical integrations is never worse than 3%. 
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Fig. 5. Modified Coulomb potential energy at B = 10 17 C Analytical formula (approximate 
P) in blue-dashed versus numerical integration (exact P) in red (the two curves practically 

coincide). 
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Fig. 6. Relative accuracy of the analytical formula for the modified Coulomb potential 

energy at B = 10 G. 

In the next section we will use the analytical expression for the electric potential energy 
(|T8l) to determine the energies of a hydrogen atom in superstrong magnetic fields. 
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4 Energies of hydrogen atomic levels originating from 
LLL in superstrong magnetic fields 

4.1 General considerations 

We are interested in the spectrum of a hydrogen atom in a superstrong magnetic field 
B. In the absence of magnetic field the spatial size of the wave function of the ground 
state atomic electron is characterized by the Bohr radius as = l/(m e e 2 ), its energy equals 
Eq = — m e e A /2 = —Ry, where Ry is the Rydberg constant. The transverse (with respect to 
B) size of the ground state of the electron wave function in an external magnetic field B is 
characterized by the Landau radius a H = 1/ \feB. The Larmour frequency of the electron 
precession is ul = eB/m e . For a magnetic field B a = e z m 2 e = 2.35 ■ 10 9 G called "atomic 
magnetic field", these sizes and energies are close to each other: as = an, E ~ u>l- We 
wish to study the spectrum of the hydrogen atom in magnetic fields much larger than B a . 
In this case the motion of the electron is mainly controlled by the magnetic field: it makes 
many oscillations in this field before it makes one in the Coulomb field of the nucleus. This 
is the condition for applicability of the adiabatic approximation, used for this problem for 
the first time in [12]. 

The spectrum of a Dirac electron in a pure magnetic field is well known [13J; it admits a 
continuum of energy levels due to the free motion along the field: 

e 2 n = ml + p 2 z + (2n + 1 + a z )eB , (21) 

where n = 0, 1, 2, a z = ±1 is the spin projection of the electron on z axis multiplied by 
two. For magnetic fields larger than B cr = m 2 /e, the electrons are relativistic with only one 
exception: electrons belonging to the lowest Landau level (LLL, n = 0, o z = —1) can be 
non-relativist ic. 

In what follows we will study the spectrum of electrons from LLL in the Coulomb field 
of the proton modified by the superstrong B . The solution, in cylindrical coordinates (p, z) , 
of the Schrodinger equation for an electron in a magnetic field B which is homogeneous and 
constant in time, in the gauge in which A = ~ B x r can be found in [TJ] . The electron 
energies are: 

\m\ + m + 1 + a z \ eB p 2 z 



E ^ m ^ -[ n P+ 2 ) m e + 2m e ' (22) 

where n p = 0, 1, 2, ... is the number of nodal surfaces, m = 0, ±1, ±2, ... is the electron orbital 
momentum projection on the z axis (direction of the magnetic field) and a z — ±1. According 
to [H], the LLL wave functions are: 

R 0m (P)= [7T(2a^) 1+ H(| m |!)]- 1/2 p l™l e (^V/(4al f )) ? (23) 

p=\p\, J \R 0m {p)\ 2 d 2 p = 1 , m = 0, -1, -2, ... 

We should now take into account the electric potential of the atomic nucleus located at 
p = z = 0. For an <C as the adiabatic approximation can be used and the wave function 
can be looked for in the following form: 

^nOm(-l) = Rom{P)Xn{z) , (24) 
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where Xn{z) is the solution of a Schrodinger equation for an electron motion along the 
direction of the magnetic field: 

1 d 2 

"2^^2 + U eff( Z ) Xn(z) = E nX n(z) • (25) 

Without screening the effective potential is given by the following formula [I]: 



Ueff(z) 

which becomes the Coulomb potential for \z\ 3> clh 

Ueff(z) 

and is regular at z = 



Rom(p)\' 



^/^Tz 2 



d 2 p , 



z~> an 



e 

\z\ 



(26) 



(27) 



U eff (0) 



(28) 



To take screening into account we must use (jTBl to modify (|26l) (see below). Since 
U e ff(z) = U e ff(—z), the wave functions are odd or even under reflection z — > —z; the 
ground states (for m = 0,-1,-2,...) are described by even wave functions. 

In Fig. 7 the different scales important in the consideration of the hydrogen atom in strong 
magnetic field are shown. 
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4.2 The shallow- well approximation 

This approximation was used in p] for the calculation of the ground state energy of the 
hydrogen atom in a superstrong magnetic field, taking screening into account. The authors 
of PP follow [5J, where this approximation was used for calculating the ground state energy 
of hydrogen atom without screening. As pointed out in [3], the accuracy of the shallow- 
well approximation for the problem under consideration is very poor. Without taking the 
screening into account it can be estimated as: 

2 



/ 



8_E_ 
~E 



\ 



ln(ln 2 H) 



1 



- 1 



(29) 



In if ) 



where if is a magnetic field measured in units of the atomic field B a , if = B/B a = (aB/aii) 2 - 
Even for B = B cr = m 2 Je = 4.41 • 10 13 G (if » 2 • 10 4 ) we obtain 5E/E « 3.5, while for 
B = 10 17 G (H = 5 • 10 7 ) 5E/E « 1.4. 

However the formula for the ground state energy in a shallow-well approximation taken 
from [5J and used in [TJ: 

~ is 



rpsw 



—2m P 



U{z)dz 



(30) 



allows us to check with which accuracy the integral of our expression ( 1T81) reproduces the 
result obtained in [T], where the potential A (z) found numerically was integrated. According 
to eq. (14) from P.R.L. paper [TJ 

-4keV , (31) 



rpsw 



-2m e a x 73. 



where E\ im is the limiting value of the ground state energy of a hydrogen atom at B — > oo. 
Substituting flU]) in (J30D we get: 
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and in the limit H 3> 37r/e 6 
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-2m P a 2 x 72.3 



(33) 



which coincides with f l3T|) within 2% accuracy (e 2 = a is used). It means that integrals over 
potential energy differ by 1% only. 

Let us present the derivation of the formula for the ground state energy of a particle with 
mass fi undergoing a one-dimensional motion in a shallow- well potential [15]. One starts 
from the Schrodinger equation: 
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Neglecting Eq in comparison with U and integrating (134"]) we get: 



X'{a) = 2p U{x)x{x)dx 



(35) 



where we assume that U(x) = U(—x) (this is why xi x ) is an even function). The next 
assumptions are: 1. the finite range of the potential energy: U(x) 7^ for a > x > —a; 2. 

X undergoes very small variations inside the well. Since outside the well x( x ) ~ e~v^^ol x ) 
we readily obtain: 

\E \ = 2fx I U(x)dx . (36) 



For 

p\U\a 2 < 1 (37) 

(condition for the potential to form a shallow well) we get that, indeed, \Eq\ <C \U\ and that 
the variation of x inside the well is small, Ax/x ~ p\U\a 2 <C 1. 

Concerning the one- dimensional Coulomb potential, it satisfies the condition ( 13 7p only 
for a l/(m e e 2 ) = a^, so (|3"6"|) cannot be used. This explains why the accuracy of (|3TJ|) is 
very poor. 



4.3 The Karnakov— Popov equation |4] 

It provides a several percent accuracy for the ground state energy for H > 10 3 . It will 
be derived in this subsection and, in the next one, we will obtain the modification of this 
equation due to screening. The Karnakov-Popov (KP) equation determines the energies 
of the families of even states with given n p and m. The starting point of its derivation is 
equation ( 125]) . Karnakov and Popov integrated it over z, from z = till a value of z such that 
an «C z as, where the assumptions of a shallow well potential (\U\ ^> \E Q \, m e \U\a 2 <^ 1) 
can be used. Assuming x( z ) ~ 1 we get: 



X\z) 



-2m e e 



—2m P e 




\Ro m (p)\d 2 p 



(3f 



Since R 0m (p) decreases exponentially for p > a H and since z is much larger than an, ex- 
panding the logarithm in powers of p/z, we can neglect the terms suppressed as (p/z) n : 
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in which 



.4 



\m\o = —m e e 



oo 



dx 



m\ 



(41) 



and the expression (1231) for Ro m (p) was used. To calculate the last integral the following 
equality must be used: 



\n(x)x ]ml e- x dx = T'(\m\ + 1) = r(|m| + l)ip(\m\ + 1) 



(42) 



where ip is the logarithmic derivative of the gamma function. 
Finally we reproduce eq. (2.11) from [3]: 

A Ho = -m e e 2 [ln(2)-^(l + H)] , (43) 

where = —7 = —0.5772, 7 is Euler constant and ip{n + 1) = —7 + Yuk=i0-/k) f° r 

n = 1,2,3,.... 

At z 3> a B the solution of the Schrodinger equation with the Coulomb potential which 
exponentially decreases at z ^> a B is given by the Whittaker function [4]. 
For the state with energy 

E = -(m e e 4 /2)A 2 (44) 
the logarithmic derivative of Whittaker function at z <C a B is [3]: 

^_ = - me e 2 S^2ln(z/a B ) + A + 21nA + 2^ - j) + ^ + 21n2 + o[{z/a B ) lnA^ J . 

(45) 

From the condition W'/W = x' we obtain an equation for the energies of the even bound 
states of the hydrogen atom originating from LLL (n p — 0, m — 0, —1, —2, ...): 



2 In — + \n2-^{l + \m\) + 0(a H /z) = 
\a H J 

= 21n + A + 21nA + 2^1--j-^ + 4 7 + 2 In 2 + 0(z/a B ) 



(46) 



and for a B 3> a B , choosing a B ^ z ^ a H , we may neglect non-calculated terms. This leads 
to the KP equation ( [4J, (2.11, 2.13)): 



ln(#) = A + 21nA + 



ln2 + 47 + ^(1 + |m|) 



(47) 



where we recall that H = B/B a = B/(e 3 m 2 e ). We postpone the analysis of this eigenvalue 
equation till the next subsection, where its modification due to screening will be established. 
The energies of the odd states are [I] : 



E, 



odd 



m e e 4 ( m^e 3 



n = 1,2 



(4f 



So, for superstrong magnetic fields B ~ m 2 e /e 3 the deviations of odd states from the Balmer 
series are negligible. 
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4.4 Equation for the energies of even states with screening 



We must equate the logarithmic derivative of wave function obtained by integrating the 
Schrodinger equation at small distances with the logarithmic derivative of the Whittaker 
function. To do this we should choose z to be larger than the distance at which the Coulomb 
potential is screened. That is why the condition z 3> an which was used in the case of the 
absence of screening is not enough. For zq > 1.5/m e the corrections coming from screening in 
(fT8"|) are only a few percent. The wave function at z > z Q is thus described by the Whittaker 
function with good accuracy. 

To derive an analog of ( 14 7p . it is necessary that x stays practically constant from z = 
up to the domain where the effective potential becomes Coulombian with high accuracy, and, 
accordingly, where the Whittaker function provides a very good approximation to the wave 
function. The characteristic distance along which the ground state wave function varies in 
the z direction is, according to the uncertainty relation: 



^char 



'2m e (-E) 



and it decreases when B increases since A ~ In 5 
with the account of it A goes to a finite value, Aoo 
of a ground state wave function we get: 



ob/A 



(49) 



'see eq. (T4T1) ). However, screening helps: 
a 11 at B — > oo (see below). For the size 



137 12 



llm P m P 



(50) 



and for 

1.5/m e < zq < 3/m e (51) 

the potential ( TT8"j) is Coulombian with good accuracy while the variation of the wave function 
[x(0) - xOo)]/x(0) is small. 

When screening is taken into account an expression for effective potential ( 126|) transforms 
into 



Ueff(z) 



\R 0m (p)f 
Vp 2 + z ' 



d 2 P 



l-e 



_|_ e - y/{2/^B+& m l z 



see Fig. 8. On the same figure a simplified potential given by 

1 



U simpiyZj 



1 _ e - V 6m i z + e 



v /(2/7r)e 3 _B+6m.2 z 



(52) 



(53) 



is shown. 
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m,,e 



-100 



DJO 




Fig. 8. Effective potential with screening for m = (blue) and m = — 1 (black- dashed), 
(Elp; simplified potential (red-dotted). The curves correspond to B = 3- 10 17 C 
Coulomb potential (green-pale) is also shown. 



The screening changes (1381) into 



= — 2m e e 2 J dx 
o 

_|_ e - v /(2/ 7r )e 3 S+6m2 a; 



1-e" 



' 6m^ a; 



|^Om(p)P 
\/p 2 + 2? 

rf 2 P = X'(^o) + Ax'(^o) 



+ 



(54) 



Calculating integrals with exponents we get: 

" Zo dx ^_ ax ^ p a dy 

Jo 



o a/p 2 + x 2 



\/y 2 + p 



2 a 2 



o \Jy 2 + p 2 a? 



In ( — ) + const 



(55) 



where, in the first approximate equality, we took into account that the integrand is expo- 
nentially damped for y > 1, and z ■>/ 6m 2 > 4; in the second approximate equality we took 
into account that pa m aijy6m| + (2/7r)e 3 £> <C 1 (because i? m is exponentially damped 
for p > an and, thus, only the values p ~ an are important). 
This leads to: 



Ax'(zo) = -2m e e 2 ln 



6m 2 



(2/7r)e 3 S + 6m 2 



(56) 



Equation fflBl takes into account the modification of the Coulomb potential at zero 
transverse direction, p = 0. In ( ]5"4"]) the modified potential at p 7^ should be used. We 



changed the factor l/\z\ occurring in eq. (JTSJ) into 1/yp 2 + x 2 in ( 154"]) . However, according 
to numerical calculations of the electric potential in the direction transverse to that of the 
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magnetic field (at z — 0) made in pQ, the screening occurs at p ^ as well. To estimate the 
role of the " transverse screening", we substituted \J p 2 + x 2 for x inside the exponents of 
(|54|) . The corrections arising in eq. (155)) turn out to be of the order of pa < a/ e 3 B/eB = e. 
The resulting correction to the ground state energy is negligible; the corrections to the 
energies of the excited states are small. In Fig. 9 we compare the effective potential given 
by fl52|) with the one obtained by substituting in the exponents of fl52l) a/ p 2 + z 2 for z for 
m = and m = —1. 




Fig. 9. Effective potentials according to ( E>|| ) i/ie replacement of z by ^ p 2 -\- z 2 in the 
exponents for m = (brown- dotted) and m = — 1 (green (pale) -dashed). For the comparison 
curves which directly correspond to 15^) are a/so shown; m = (blue) and m = — 1 
(black- dashed). The curves correspond to B = 3- 10 17 G. 



With the help of f )56|) we get the modified KP equation, which takes screening into 
account: 



In 



H 



\ 



1 + h H J 



A + 2 1 1 ) A + 2 ( ' ( 1 — - ) In 2 \ h 1 rf] + |/;/|) . A7) 



The ground state has m = 0; for H = 10 3 , its energy as given by ( T57|) coincides with the 
result of the numerical solution of the Schrodinger equation at a few percent level. 

For H > 10 5 the expansion ip{l — i) = -0(1) — |^ has better than 1% accuracy [1]. Using 
this approximation, we obtain: 

4 

H = 10 5 ; 5 = 2.3 • 10 14 G ; A = 6.91 ; £ = ^— A 2 = -0.65 keV . (58) 

The accuracy of the adiabatic approximation is characterized by the ratio oh/ob = ^/B a /B 
which equals 0.0032 for H = 10 5 . 

The screening starts to be important at H ~ 37r/e 6 ~ 10 7 . For example, H = 10 7 and 
H = 1.7 • 10 7 correspond to the same ground state energy respectively without and with 
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screening. For H ^> 10' 
the limiting value 



which corresponds to B ^> 10 16 G, the ground state energy reaches 



A 



lim 



« 11.2 , M im « -1.7 keV , (59) 

which differs from the result obtained in [1] by a factor 2.5. So, when screening is accounted 
for, the accuracy of the shallow-well approximation is never better than 250%. 

For each value of m (m = 0, — 1, — 2, ...) eq. ( 1571) gives a tower of even states, starting 
from the ground state, the energy of which is very large on the Rydberg scale. Let us rewrite 
(1571) for the m = tower: 



/ 



In 



H 



1 + —H 

3tt 



A + 21nA + 2V>(l 



In 2 + 37 



(60) 



The ground state energy corresponds to a large value of A, of the order of the left hand side 
of ( 160]) as just discussed (see ( |58l) . (1591)). ip(x) has poles at x = 0, —1, —2, ... which correspond 
to A = 1, 1/2, 1/3, ... The Balmer series of energies E n = — m e e 4 /(2n 2 ), n = 1, 2, ... that they 
form would correspond to an infinite value for the left hand side of (l60j) . That it is always 
finite shows that even states get shifted from their Balmer values. 

Let us determine the energy of the first excited even state. From the well-known ex- 
pansion of the gamma function at small values of its argument we obtain the corresponding 
expansion for ip: 

T(e) = -- 1 + 0(e), V>(-e) = --7 • (61) 



Plugging in (160|) A = 1 — e we obtain: 



A = 1 - 



(62) 



In 



H 



\ 1 + * H ) 



- In 2 - 7 + 1 



Eq. ( )62|) describes how the first excited even state (the second even state) moves towards its 
Balmer value A = 1 with growing H. When H — > oo we obtain A 2 ^ = 0.78 which determines 
the energy of the first excited even state in the limit of infinitely large B according to the 
formula E = 



m, 



>e 4 /2)A 2 . The first odd state has the energy 



odd 



-m e e 



72- 



For arbitrary n the analog of (162]) looks like 



A, 



2/n 2 



n 



( 



(63) 



In 



H 



- In 2 - 7 + 1/n + 2 In n - 2 YJlll 



n—l 1 
k 



where the last sum in the denominator is the harmonic number i? n _i. 

A similar formula with no screening and in the limit H ^> 1 can be found in [16] ; we 
however disagree concerning the factor 2 in the numerator of their eq. (11). An expression 
which coincides with the leading term in ( 163]) for H 3> 1 without screening is contained in 
the book by Khriplovich "Teoreticheskii kaleidoskop" , Novosibirsk, 2007, eq. (3.39). 
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The energy spectrum of the hydrogen atom in the limit of infinite magnetic field is shown 
in Fig. 10. 
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Fig. 10. Spectrum of hydrogen levels in the limit of infinite magnetic field. Energies are 

given in rydberg units, Ry = 13.6 eV. 



5 Conclusion 

An analytical formula for the Coulomb potential &(z) in a superstrong magnetic field has 
been derived. It reproduces the results of the numerical calculations made in [1] with good 
accuracy. Using it, an algebraic formula for the energy spectrum of the levels of a hydrogen 
atom originating from the lowest Landau level in a superstrong B has been obtained. The 
energies start to deviate from those obtained without taking the screening of the Coulomb 
potential into account [I] at B > 3irml/e 3 ~ 6 • 10 16 G which is much larger than the 
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largest magnetic field known to exist today, that of neutron stars called magnetars (where 
B m 2- 10 15 G may exist). Without the presented explicit calculations we would not know at 
which B substantial deviations of the energies of atomic levels from those obtained without 
taking screening into account take place. 
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Appendix 

Electron mass in a strong magnetic field 

In the non-relativistic approximation, from (|2~T]) we obtain: 

p 2 , I.eB eB , . , . 

e n = m e + -£*- + (n+-)— + a z —- , Al 
2m e 2 m e 2m e 

where the last term originates from the interaction of the electron magnetic moment with 
the magnetic field. Taking into account the anomalous magnetic moment of electron for the 
energy of the lowest Landau level (n — 0, a z — — 1) we get: 

p 2 a eB , . , 

s = m e + ^-- . A2 

2m Air m e 

This result is valid for a magnetic field B < m^/e, while for B S> m 2 e /e the correction to 
the energy changes sign and the strong power dependence on B is replaced by a double 
logarithmic one |17j : 

2 

e = m e + 7 ^ + -^m e \n 2 (eB/m 2 e ) . (A3) 
2m P Air 



As a result, for the energy of the ground state of the hydrogen atom, we obtain: 



E Q = m e 



1 - —A + — In 2 — 

2 47r ml 



(A4) 



At H pa 10 5 (see ( 158]) ) the last term in brackets compensates the second term and Eq 
becomes bigger than m e for a bigger H . However in the differences of the energies of atomic 
states originating from the LLL the universal correction (the last term in brackets in (1A4|) ) 
cance lfl and for the "experimentally observable" energies of transitions between these states, 
the formulas obtained in Section 4 should be used. 

Concerning the double logarithmic term in flA3|) it was noted in [TB] that the appearance 
of an "effective photon mass" at fields B > m^/e 3 leads to the transformation of the double 
logarithmic dependence on B into a single logarithmic dependence (concerning extraction 
and summation of higher order terms in a \aB see also [19], [20J). 



7 We are grateful to V.S. Popov for this remark. 
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